function printout_dispersion_equation (m, CONSTS, plot_data)

    p_high = 25; %15
    num_dots = 2000;
    p_vec = (1.0 : p_high/num_dots : p_high)';
    
    for variant = 2 % 1:2
    
        if (variant == 1)
            y_vec = zeros(size(p_vec));
            for i = 1:size(p_vec, 1)
                p = p_vec(i);
                y1 = matrix_of_boundary_conditions(p, m, CONSTS, false);
                y_vec (i) = det(y1);
                %fprintf('p=%d    %e + %ei \n', p, ...
                %                 real(y_vec), imag(y_vec));
            end
        elseif (variant == 2)
            y_vec = func_delta(p_vec, m, CONSTS);
        else
            error('Wrong variant');
        end

        if (plot_data)
            figure; plot(p_vec, real(y_vec), 'b.-'); title('real(determinant)');
                     xlabel('p'); ylabel('\Delta');
            figure; plot(p_vec, imag(y_vec), 'b.-'); title('imag(determinant)');
                     xlabel('p'); ylabel('\Delta');
        end

    end
end